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Abstract: We study two constrained scalar models. While there seems to be equiv- 
alence when the partially integrated Feynman path integral is expanded graphically, the 
dynamical behaviours of the two models are different when quantization is done using 
Dirac constraint analysis. 



INTRODUCTION 



There are many models in the literature which claim that they are equivalent to Quantum 
Chromo dynamics (QCD) in some sense. We can site a recent work by Hasenfratz and 
Hasenfratz Z 1 , or the work of Akdeniz et al / 2 among many papers concerning this topic. 
We mention these two papers since in these papers the starting lagrangians look very 
different, although the effective lagrangians obtained after some manipulations are made, 
are exactly the same. There are many other papers in the same spirit. One example is 
reference 3, which , in some sense, gave rise to reference 1. Another example is reference 
4, which was actually a pioneering paper in this endeavor of finding models whose effective 
lagrangian looks like the standard model. 

It was always a puzzle to us how similar manipulations made on very differently 
looking Lagrange functions resulted in completely the same effective lagrangian. In this 
note we try to investigate this phenomena using scalar models. We think that our results 
in the scalar case may give additional information on this phenomenum. 

We will use two constrained scalar models. In reference 1, the authors imposed the 
constraint = ^i^^T a il) = on the fields of the free spinor lagrangian. To resemble 
this model here we first study the case where we write a lagrangian which is essentially 
equivalent to L = ^d^d^cp . We will impose a different constraint, though, since a 
constraint </><9 M = results in a truly trivial model. We can also introduce inner symmetry 
to the theory and make a O(N) model along similar lines. For the time being we do not 
pursue this . 

The authors in reference 2 impose the condition that their current equals product 
of vector fields instead of zero, = A^A 2 , which differs from the constraint used in 
reference 1. This complicates the problem, but all of additional fields introduced to the 
model decouple and at the end only one vector field survives. The propagator for this field, 
and only for this field, is generated in the one loop calculation. At this point the resulting 
effective theory looks exactly like that of reference 1. 
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We have doubts whether these two models are actually equivalent to QCD in all 
aspects. One may refer to an old work of Wilson / 5 and to a more recent work of Zinn- 
Justin / 6 , and using the calculations made in reference 7, claim that these two models are 
actually examples of trivial models / 8 . 

We will not dwell on these points here. We will only investigate in what sense two 
models are equivalent when the effective lagrangians derived from them seem so. In the 
next section we present two constrained scalar models. We get a theory which is totally 
trivial if we impose the current made out of scalar fields equal to zero, the analogous case 
as given in reference 1. We instead use two models where the current is equal to one and 
two auxiliary fields, thus introducing eight and sixteen new degrees of freedom respectively 
plus constraints that will eliminate these . We study the Dirac bracket relations satisfied 
by the respective fields . We see that the new introduced vector fields via the constraint 
equations somehow replace the canonical momentum of the scalar field. 

In Section III we derive effective lagrangians for these two cases and show why do 
they seem to be equivalent on this level. We end with some remarks. 

II. Quantization of the Models using Dirac Constraint Analysis 
II. A 

We start with 

L A = l -d^4> + igX^d^ + \g 2 <f\ 2 . 

We are in four dimensional Minkowski space and \i takes the values zero to three 
is an auxiliary field with no kinetic term, g is a coupling constant. 

The equations of motion are 

«W + igd^<f>) = ig\»d»<j> + g 2 \ 2 ^ 2 

igcpd^cp = -g 2 \^4> 2 3 
2 
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which can be shown to be equivalent to 



d^(p = 0. 4 

In this calculation we use the methods given in Dirac's book / 9 . The canonical mo- 
menta are 

ti> = (do + ig\o)(j) 5a 

and 

vta m =0 5b 

which gives us four primary constraints. 
The canonical hamiltonian is 

H=-(ir <f> - igX (P) 2 + igX^dicj) + -d^d^ - ^ 2 A 2 2 . 6 

We get secondary constraints when we set the Poisson bracket of H with 7Ta m equal to zero. 

= ^tta m = (igfat + g 2 Xo4> 2 )go ll - ig^d^g^ + g 2 X^<p 2 =0. 7 

We take 

H E = H + c m 7Ta m 8 

and further calculate 

j t K^ = [H E ,K^ 9 

where square brackets mean Poisson brackets. Note that Ao appears in Kq, and the A, in 
Ki. We get one equation with Co when [He, Kq] is calculated, which fixes the value of cq 
and does not give additional constraints on the system. [H E ,Ki] give us equations which 
fix Q . We do not get any additional constraints. 

We can calculate the Poisson brackets between the different constraints. 



[tta , Kq] = 2g 2 (p 2 g 00 , [0a* , Kj] = -g 2 (p 2 g ij7 
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[K ,Ki] = g 2 c\>d l <\>-2ig 3 X^ 2 

all the other brackets of the constraints with each other are zero. We see that all these 
brackets are second class. We calculate Dirac brackets between different fields. 

[cf ) ( x ) 1 X (y)] D = -^-5 s (x-y), 10a 
2g<p 

{<p(x),X t }=0, 10b 
which shows that A is like 7r^, and \ decouples. We can set \ equal to zero . 



II.B 

We propose another model where 

Lb = \d^4> + ig(X^ + A^d^ - gA^A 2 11 
The primary constraints are 

7T A = 12 

7T A =0 13 

The hamiltonian reads 

1 1 

H = -(ti> - ig{\ + A )4>) 2 + -d^d^ + ig(X t + A^d^ + gX^A 2 14 



where 



n<f, = d <f) + ig(X + A )(p 15 
H E = H + c^+c^ 16 



Secondary constraints 



[H E ,7r A ^] = Ql = 17 

[H El 7t x J=Ql = 18 

are given as 

Q]i = ig(f>(n<t> - ig{X + A )(p)g fl o - igipdifaifj, - gX^A 2 - 2gX v A v A^ 19 



Ql = W<t>{K<t> - ig(>>o + A )(p)g^ - ig(pdi4>gi^ - gA^A 2 20 

We see that the system is closed, since the Poisson brackets of Q^Q^ with He involve 
eight coupled equations for and c 2 . We get no further constraints. 

When we calculate the Poisson brackets of the constraints with each other, we see 
that they are all of second class. We have sixteen second class constraints and eighteen 
degrees of freedom. We have traded some of our fields in terms of others, but we did not 
change the number of independent variables. 

Now we can calculate the Dirac brackets between different fields. We are particularly 

interested in the brackets between <fi and A M , A M , since the Poisson brackets between the 

same fields are zero. The effect of the constraints in the system are reflected to the Dirac 

brackets; hence, they do not vanish when thay are taken between and the auxiliary fields. 

We give below the result of some sample calculations. 

rj , , „ , „ Da , , i<j) 2 (-2A A 1 )[g<j) 2 (A 2 -2A 2 ) - 3A 4 ] _, 3W 
[<P(x),A 1 (y)] D 4>(y) = nW \ ^\x-y), 21a 



(x),A (y)] D cf>(y) = 



^-2AW(^-2A 2 )-3^ V)(x _ m 



[m.uvrm = ^ a -^ a )[w-2A a )-^-4A^] J( , )(a _ y)t 21c 

Mx)Mv)\ D m = 2ift^M^ - AA 2 A fl M l ) - M X 2 + A 2 X o)] tf(s)(g _ y) ^ 
where 

A = [g 2 (p 4 (A 2 - 2A 2 ) 2 - (6A 2 + A\^)(A 2 - 2A 2 )A 2 g( p 2 + 9A 8 ] 22 
Here A 2 means the three vector A squared. 

Upon quantization we see that A^ and A M seem to contain part of ti^. We expect only 
Tifjj to have nonzero commutation relations with <p and in this model both A^ and A M also 
will have non zero commutations with 4>. The constraints Q 1 ^ = and Q 2 ^ = relates tv^ 
to these fields. 



The model we studied seems to be considerably different from the one studied in the 
first Section. The fields in the model has non zero Dirac brackets; so, we can not set them 
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equal to zero, as in the previous model. The space components of the vector fields do not 
decouple and can not be set to zero. 



Note that in both of these models the degrees of freedom is two. In Model A we start 
with ten degrees of freedom, two for the <j> and eight for the A field and their respective 
momenta. Eight constraints reduce these to two. In Model B we start with eighteen 
degrees of freedom since we have two vector particles. Sixteen constraint equations reduce 
this number to two. As far as the equations of motion are considered these two models do 
not seem to be alike. 

III. Feynman Rules using the Path Integral 

Here we study the two models using Feynman diagram expansions of the path integral 
after the integral is partially integrated. We start by studying model A, then contrast our 
results with that of Model B. 

III. A Here the path integral is written as 

Z = J dtpdn^dX fl dn\ ll 6('irx ll )6(K fl )detM tlu expiS, 23 



where 



24a 



S = J d 4 x[n^d (j) + tta^qA^ - H E ]. 246 

We write the Dirac delta functions in the integral form, introducing new variables and 
express the determinant in the exponential form using ghost fields. 



detM^ = J dc+dc u e* c i M ^ c ». 
The integrations over the momenta and are performed easily and we end up with 

Z = N J d^dA^c+dc^-*'^-^^^ 

25 



where we define iV M = A M — A^. We can calculate the inverse propagator, D~l for the N p 
field by taking two derivatives of eq. (25) with respect to the iV M field. In the momentum 
representation we get 



n -i (n \ = _ n 2 f d4 P ( Pn + <ln)(Pv + 2 qv) 
^ yq) y J {2%Y K 



r) 4 p 2 (p + q) 2 

which looks like the massless vector boson propagator, at least in a particular gauge. Note 
that all the components of the vector field have non-zero propagation . 

All other fields have zero propagators if we use dimensional regularization. Here we 
set J d^p^s = 0. When we drop all the fields with zero propagators we end up with 

Seff = ~Trlog(-d 2 + igN^ - igd p N^). 27 

Upon expanding the logarithm we can evaluate the multi-point functions for the N p fields. 
Eq. 26 dictates a necessary condition on the coupling constant g , though, to have a well 
defined expression for the propagator function given by this equation, which reads 



This condition makes the model asymptotically free in the ultraviolet regime. 



By taking all the non vanishing terms we see that for the composite field A M the effective 
lagrangian can be written as 

L eff = \d^N v d^N u + d^N„d u N» + gf^ p NPN v N p + g 2 V^ pa N»N v N P N° . 29 

Here f^ vp is proportional to momentum and Kronecker deltas and V pvpa is made out of 
Kronecker deltas. Higher order functions, starting with the fifth point function, drop with 
higher powers of g. For example the five point function goes as g 5 . They do not fit into 
this scheme of effective lagrangian and are calculated as loop corrections. 



Here we calculated the Feynman rules for this model and showed that apart from the 
restriction dictated by eq. 28, we get rules similar to those as a gauge theory. One can 
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calculate physical processes using these rules and will get free parton model results, as is 
the case in a similar model I 1 due to the restriction dictated by eq. 28. All the physical 
processes that involve interactions will involve powers of the coupling constant which goes 
to zero . Any possible divergences due to loops will be canceled by the zeroes coming from 
extra powers of the coupling constant. Only terms which do not involve any interactions 
are finite. These terms are the same as those given in the free field case. 

III.B The path integral for Model B, in the hamiltonian formalism, is written as 

J dA^dn A ^dnx^dX^d(pdn ( f ) d(nx l ,)S(7VA fl )S(Ql)S(Ql)(detM)expiS 30 

Here 

rf 4 x[7r^5 (/) + 7Ta m <9 A m + tta^oA^ - - + ig(X + A )(f)] 2 - -d^d^ 

-igiK + AA^-gX^A 2 ] 31 
Ql = 0(tt - (A + A o )0)^o - Wifoin - A M A 2 - 2X v A v A li 32 
Q\ = <f>(^4> ~ ( A o + A )(p)g fl o - 4>di(pgi^ - A^A 2 33 

M is a eight by eight matrix whose entities are made out of derivatives of and Q 2 ^ with 
respect to the fields A M and A M . 

We can use the integral representation of the Dirac delta functions. 

KQl) = ^J dB^exp-iB^Ql 34 

KQl) = ^J dE^exp-iEfQl 35 

Using ghost , i.e. Grassmann valued fields c M , e M , c^, e^, we can raise detM to the expo- 
nential. 

detM = J dc^dcvde^depexpiN 36 

where 

N = (<i + ej t )(^ 2 2 ^o^o)(c i , + e„) + cl^gA^K - 2gg^X K A K - 2gX^A v - 2gX„A fl )c u 



+ct (-g^A 2 - 2gA»A v )e v + e^-gg^A 2 - 2gA»A»)c v 37 
When the momentum integrals are performed we get 

Leff = i[\d^4> + igG^cj) - 9 - (B + £ ) V - 9^ A 2 + gB^X^A 2 

+2gB»A*\ v A v + gE^A 2 + g 2 f^ 2 f + c^-gg^A 2 - 2gA»A»)f v 
+fl{-gg^A 2 -2gA»A v )e v +c\[2gA 2 g^-2g\»A v -2gg» v ^ 38 
Here / M = + e M . We set = A M + - B^ - E^. 

We perform the integration over and obtain 

S eff = -\Trlog[-d 2 + igG^ - igd^G -^(B + E ) 2 + ^/oVo] 



+ J d 4 x [—gX^A^A 2 + gB^A 2 + 2g\» A^A V B V + gE^A^A 2 + e\(-gA 2 g^ - 2gA' t A v )f v ^ 

+fl(-gA 2 sr - 2gA"A v )e v + c\(2gA 2 g^ + AgA»A» - 2g\» A" - 2g\ v A» - 2g\ p APg^\\ 

39 

Note that only G M propagates among all the fields given above. To find the propagator 
we take two derivatives with respect to the respective fields. 

d 2 S eff -1 , l 4 (p» + q »)(p» + 2q») 

lo = ~ \j9 / a p ; .„ 40 



dG^x)dG v (yY" (2vr)^ J " p 2 (p+q) 2 

Subscript zero on the derivative means that all the fields are put to zero after the differ- 
entiation is performed. 

Note that this is the same expression for the propagator of the A M field as given in 
eq.(26) We also see that 



d 2 S e ff _ d 2 S e ff _ d 2 S e ff 

~dBf ~ dE 2 ~ d gldg Q ~ (2 



i r d^p 



41 



This expression is zero by dimensional regularization.All the other fields also have zero 
propagators since the effective lagrangian does not have any terms which are only bilinear 
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in these fields. All these terms involve quartic interactions of these fields. When we drop 
all the field with zero propagators we end up with 

S'Jff = -^(-d 2 - igG^ + igd»G») 42 

This is the same expression we found for Model A. Therefore all the results obtained for 
Model A from this expression are also true for Model B. We can not differentiate Model 
A from Model B as far as perturbative expansion in terms of Feynman diagrams are 
concerned. 



Conclusion 

Here we have studied two very dissimilar models which have the same Feynman expansions. 
A complete constrained Hamiltonian analysis shows that the two models are different. 
One reason we have studied this problem is to be able to clarify the behaviour of many 
Nambu-Jona-Lasinio/ 10 ' 3 like models which are claimed to be similar to QCD Z 11 . There 
are people Z 12 who disagree with this equivalence. The claim in reference 12 is that after 
an investigation of a lattice Nambu-Jona-Lasinio model both by the Monte Carlo method 
and Schwinger-Dyson equations, studying renormalization group flows in the neighborhood 
of the critical coupling where the chiral symmetry breaking phase transition takes place, 
in no region of the bare parameter space renormalizability of the model is found. We 
propose that, in addition to the standard methods of looking at the renormalization flow 
and fixed point structure of two models to show equivalence, their constrained analysis 
may be another check. Still another method is to study the predictions of these models for 
different physical processes . An old calculation / 7 and old paper / 5 seem to suggest that 
the Nambu-Jona-Lasinio type models may indeed be trivial at four dimensional space-time, 
perhaps like the 4 model is. 
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